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We study the number of propagating Bloch modes Nb of an infinite periodic billiard chain. The 
asymptotic semiclassical behavior of this quantity depends on the phase-space dynamics of the 
unit cell, growing linearly with the wavenumber k in systems with a non-null measure of balHstic 
trajectories and going like ~ in diffusive systems. We have calculated numerically Nb for 
a waveguide with cosine-shaped walls exhibiting strongly diffusive dynamics. The semiclassical 
prediction for diffusive systems is verified to good accuracy and a connection between this result 
and the universality of the parametric variation of energy levels is presented. 
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I. INTRODUCTION 

According to the Bohigas-Giannoni-Schmit conjecture 
[T], a system with chaotic phase-space dynamics is 
expected to show universal quantum level statistics in 
the semiclassical limit, consistent with the predictions 
of Random Matrix Theory (RMT). For instance, the 
two-point level correlation of closed chaotic systems has 
been shown experimentally and numerically to agree 
with the Gaussian Orthogonal Ensemble (Gaussian Uni- 
tary Ensemble) for systems with (without) time-reversal 
symmetry |21 |2j- Analytical results supporting this 
universality have also been achieved using the semi- 
classical trace formula, first by means of the diagonal 
approximation and more recently including non-diagonal 
terms verifying full agreement of the semiclassical form 
factor with RMT g]. 

Much less is known for spatially extended chaotic 
systems. In the case of systems that relax by a diffusion 
process, the spectral properties deviate from random 
matrix theory in a particular way. Dittrich et al. [SI |S] 
considered the statistical properties of the energy bands 
Enfi of a closed periodic diffusive system with N chaotic 
unit cells and analyzed the signatures of chaotic diffusion 
in the form factor (the Fourier transform of the two-point 
level correlation function). They showed that their semi- 
classical results agree nicely with those of Simons and 
Altshuler [7i |8| who proposed that the fluctuations of the 
energy bands E^^ or in fact any parametric variation 
of the energy levels are universal within the appropriate 
Dyson ensemble (depending on time-reversal symmetry) . 
Simons and Altshuler calculated the universal two-point 
level correlation function for a disordered system under 



the influence of an external perturbation -for instance 
in the form of an Aharonov-Bohm flux- and showed 
that the response of the spectrum to those parameter 
variations was universal after a suitable scaling with 
two system specific parameters. Dittrich et al. [SJ |B] 
identified one of these parameters with the diffusion 
coefficient of the chain. 

These correlation functions or form factors provide 
a detailed description of the statistical properties of 
the spectrum. In this work we consider a particular 
quantity directly related to the correlation function 
and of great practical importance for wave propagation 
in periodic media: the number of propagating modes 
Nb{E) with energy E. In a periodic system the prop- 
agating modes are given by the Floquet-Bloch modes, 
which are quasi-periodic solutions of the Schrodinger 
equation. They propagate ballistically through the 
system with constant (group) velocity and are the 
only modes allowed to transport energy. For instance, 
the dimensionless Landauer's conductance of a long 
periodic sample fluctuates closely to the number of these 
modes propagating in a given direction. In Faure 
studied this quantity for quasi-one-dimensional periodic 
systems in the small Plank constant limit h — 2TTh — 0. 
He showed that if the phase-space displays tori in 
the transverse stroboscopic Poincare section, corre- 
sponding to trajectories that propagate ballistically 
through the system, then Nb ^ h^^ but, on the other 
hand, for fully chaotic and diffusive systems Nb ^ /i^^^^. 

We study Nb{E) numerically for periodic waveguides 
with diffusive classical dynamics and finite horizon in the 
semiclassical regime. The paper is organized as follows. 
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In Section |TT] we review Faure's derivation of the number 
of propagating Bloch modes in a waveguide. In section 
|III| we define the bilUard where the predictions are tested. 
The results are summarized in section llVl Then in section 
W\ we discuss in further detail our results and its relation 
to previous works and present some conclusions. 

II. NUMBER OF BLOCH MODES FOR fi ^ 




In order to be self-contained, we start reviewing the 
essential results of [^. First, in an infinite periodic 
waveguide, the distribution of quantum velocities 
Pfi{v,E) determines the number of propagating modes 
Nb{E) with positive velocity [see Eq.S]. Second, the 
asymptotic form of the average {NBiEy) for /i — ^ 0, ob- 
tained using a semiclassical approximation of {Pri{v,E)) 
[see Eq. 



FIG. 1: Part of the spectrum k„{9) for the bilhard studied in 
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Consider a periodic waveguide in the cartesian (x, y) 
plane, with x g] — oo, +oo[ (the propagation axis) and y 
bounded for every x between the wave boundaries that we 
define as hard walls. The waveguide is composed of cells 
with length L. In the waveguide particles are free, hence 
Schrodinger's equation reduces to Helmholtz equation, 
i.e. 



(V2 + fc2)^(a:,y) = 



(1) 



with Tpi^j y) satisfying Dirichlet boundary conditions at 
the walls. The reader can have in mind several physical 
contexts where this is valid such as optical, microwave 
or acoustic cavities or degenerate electron gases in meso- 
scopic devices. The classical analog of the quantum sys- 
tem under consideration is given by the free particle 
Hamiltonian H{x,Px,y,Py) — {pi +Py)/2 (we consider 
particles of unit mass), plus perfect elastic walls as the 
waveguide boundaries. The results we present below hold 
for general periodic Hamiltonian systems (not only bil- 
liards) where the confining potential in the transversal 
direction y could be soft. Eq. (13 1 is the general result 



for that case whereas the particular result for hard-wall 
billiards is given by Eq. ( 15 ) which is the expression we 



use in the rest of the paper. 

Let us consider the phase-space surface of constant en- 
ergy in one unit cell Y,e = {{x,Px,y,Py) / < x < 
L, hi{x) <y <h2{x), H{x,px,y,Py) = E} where hi{x) 
and /i2 (x) define the lower and upper walls of the waveg- 
uide respectively. Since the boundaries of the guide are 
periodic and S^; is compact, we have from Bloch the- 
orem that the solutions of Schrodinger's equation [Eq. 
([T|] form bands, i.e., its spectrum has the form E^^e with 
n e N the band index and 9 = qL G [— tt, tt] where q is 
the Floquet-Bloch wavenumber. Figure [T] illustrates this 
spectrum for a chaotic periodic waveguide. An important 
consequence of Bloch theorem is that the propagation in 
a periodic system is ballistic with a group velocity [TU] 



L dE, 

n 



de 



(2) 



Sec IV Note that for the statistics of Nb (k) the computation 
consider values of k considerably larger than the ones of this 
figure. 



Since N = jl^dx6{C - f{x))Q{f' {x))^ is the number 
of points such that f{x) — C with positive derivative in 
the interval (a, 6), Faure expressed the number of Bloch 
modes with energy E propagating with positive velocity, 
Nb{E) as 



Nb{E) = Y,I de 8{E - E^^g)Q{v,,^g 



dEr, 



de 



(3) 



where 6{-) is Dirac's delta function and 9(-) is Heaviside 
step function. Considering the distribution of quantum 
velocities Pti{v,E) defined as [51 [TT] 



Pniv,E)^h^J2 J ^5{v-Vnfi)5{E-En,^ 



with Vn_0 defined in Eq.(p|) one obtains 



Nb{E) 



1 



vPn{v,E) dv. 



The normalization over v of Pniv, E) is 



de 



dvPn{v,E)^h^ I ^piE, 
Zn 



(4) 



(5) 



(6) 



with p{E, e) = '}2y^{^ ~ En,e) the density of states. 
We note that Eq.jsl) is an exact result, a particular 
case of the Kac-Rice formula [12], valid for quasi-one- 
dimensional periodic systems [9]. 

The number of propagating Bloch modes is a rapidly 
fluctuating quantity as a function of E at scales near the 
mean level spacing AE — 1/pei where is the mean 
density of states. These fluctuations are eliminated by 
considering the average of Nb{E) over an energy inter- 
val SE such that E '3> 6E :$> AE (i.e., quantum me- 
chanically large but classically small), which is the same 
definition used to determine pE from p{E,e). This pro- 
cedure is supposed to play a similar role to the ensemble 
average for disordered systems. Note that we can write 



{Nb{E)) = 



1 
'hi 



v{Pn{v,E))dv. 



(7) 
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As follows from Eq.([6]), {Ph{v,E)) is normalized to 
hP'pE = ve where ve is the Liouville measure of the 
constant energy surface (Weyl law). 

The expression at the right hand side of Eq.Q will 
be used for the theoretical semiclassical analysis of 
{Nb{E)) but not as a recipe to implement a numerical 
computation because it is difficult to properly sample 
small velocities to generate the velocities distribution. 
However, we are not concerned with the numerical 
evaluation of this distribution but with the average of 
Nb{E) which can be computed in different ways. In 
Sec. IIII BI we describe a numerical method that allows to 
compute {Nb{E)) accurately bypassing the computation 
of Pn{v,E). [See in particular the last paragraph at the 
end of the aforementioned section] . 

We now turn to obtain the leading order term of the 
semiclassical expansion of {Nb{E)). Asch and Knauf [TT] 
proved the asymptotic convergence of {Ph{v, E)) towards 
the classical asymptotic velocities distribution Pa{v,E), 



Pniv, E) = Pa{v, E) when H^Q, 



(8) 



that holds for test functions independent of h, i.e. over 
intervals of order 5v ^ 1, SE ~ 1, a coarse graining 
consistent with Q. In the right-hand side of ([s]) Paiv, E) 
is the probability density of the asymptotic mean velocity 



lim 



x{t) 



(9) 



where af = DE/t. In the limit i — >■ oo, Pt{v, E) converges 
to a punctual measure at v = 0. It has been shown [13], 
using Gutzwiller trace formula and also with RMT argu- 
ments, that in systems with hyperbolic classical dynamics 
the variance of a quantum operator in the semiclassical 
limit is equal to the variance of the associated classical 
observable at time t = tfj/g, with = h/AE = 
the Heisenberg time of the unit cell and g a factor de- 
pending on the anti-unitary symmetries of the system. 
For the billiard we study in Sec. |III| the factor is g = 2 
due to the transversal mirror reflection symmetry of the 
unit cell [55] (see the discussion in Sec. Iv|. Then, we 
have a relation between the variance =7'^^ g) of the 
quantum level velocity and the variance of the classical 
mean velocity at 



gPE 

tH 



for -> 0. 



(11) 



This relation between level velocity variance and the 
classical diffusion coefficient De for the energy shell I^e 
was studied in detail for the kicked rotor [T3] where the 
diffusion coefficient displays a nontrivial behavior as a 
parameter is varied. In a different context |15j explores 
the relation between classical diffusion coefficient and 
quantum evolution properties. 



Inspired by Eq.(ll), Faure conjectured a higher order 



semiclassical equivalence between the averaged distribu- 
tion of quantum velocities {Pri{v,E)) and the classical 
mean velocity distribution at time t = tn/g, i.e.. 



where x{t) is the longitudinal position at time i of a free 
particle in the waveguide whose initial condition was 
taken randomly with a uniform probability distribution 
on the surface T,e- When the classical dynamics is 
not purely ergodic, so there are some tori associated 
to ballistic trajectories in the stroboscopic Poincare 
section, the leading order term in the semiclassical limit 
of Nb can be obtained just using equivalence ([8| in 
the integrand of Eq.([7]). On the other hand, in case 
the dynamics is completely ergodic, Pa{v,E)dv is a 
punctual measure at w = because Va — almost 
surely. To obtain the leading semiclassical contribution 
in {Nb{E)), it is necessary to quantify more precisely 
how the quantum level velocity distribution approaches 
the classical asymptotic velocity distribution as fi — > 0, 
that is, we need a refinement of ([s]). 

When the classical dynamics is fully ergodic and the 
mixing rate is rapid enough, we expect that (almost) any 
initial ensemble of particles in the system will relax by 
diffusion, with a Gaussian density profile spreading ac- 
cording to {x{tY) = DEt for large times, with De the 
diffusion coefficient for Yie- Hence, the mean velocity 
distribution at time t is given by 



Ptiv,E) 



VE 



2'Kat 



exp 



2af 



(10) 



Pn{v,E)^Pt^t^/g{v,E) when h 



0, 



(12) 



over widths of order SE ^1, Sv ^ Vh. This conjec- 



ture is based on result (111 and the Quantum Ergodicity 
Theory [TB] . It can also be seen as a consequence of the 
universality in the parametric variation of energy levels 
discovered by Simons and Altshuler [51, as we discuss in 
Sec. |V] Using (12) in the integral of Q we have, for a 
purely diffusive waveguide. 



{Nb{E)) = 




(13) 



Equation ( 13 ) was derived in [5] without the factor g 



and checked numerically for the kicked Harper model. 

Billiards are a particular class of systems with the sim- 
plifying property that their phase-space dynamics is the 
same on every energy shell except for a change of time 
scale. For two-dimensional billiards fE = 2TrAc depends 
only on Ac, the area of the unit cell. Hence, the semi- 
classical limit /i — is equivalent to A: — > oo due to 
the relation fi^k^ = 2E. In S^; the upper bound of the 
mean asymptotic speed is \J2E. Considering the change 
of variables v = \j2Eu we have (from Eqs. [t] and |8| 
{Nb[E)) = J^uPa{u)du with Pa{u) the classical 



FIG. 2: Portion of the infinite cosine waveguide with a tra- 
jectory's segment inside. 



asymptotic velocity distribution in S1/2 where the parti- 
cles move with constant unit speed. Since \/2E — hk, 



{NB{k)) = 



k 
2^ 



vPa{v) dv. 



(14) 



Thus, for balHstic waveguides {Pa{v) ^ VESiv)) the aver- 
age number of propagating Bloch modes goes asymptot- 
ically for k going to infinity as ~ ^baik/2TT + o{k) where 
/iba/ is the measure of ballistic trajectories in a strobo- 
scopic cross section of the waveguide as defined in [3] . For 
diffusive billiard chains, the diffusion coefficient satisfies 
De = \v\Di with speed |w| = \/2E and Di the diffusion 
coefficient in the energy shell S1/2 where particles move 
with unit speed. Then, using |ti| — hk and ve — ^nA^ in 
Eq.( 13 1 we have that the average number of propagating 



modes in a diffusive wave guide is 



2ttL 



^ + o{Vk) 



(15) 



In the following sections we will study how well this re- 
sult applies to a version of the cosine waveguide, where 
many semi-classical results for quantum systems have 
been tested flTllTSl. 



III. COSINE BILLIARD CHAIN 
A. Definition and classical dynamics 

We consider a periodic bidimensional waveguide with 
cosine-shaped hard walls [H]. The unit cell is defined 
for — 1 < 2; < 1 with the upper boundary given by 
h2{x) = 1 + ^(1 + cos(7ra;)) and the lower boundary 
defined by hi{x) = ^{1 + cos(7ra;)). The waveguide is 
constructed as an infinite one-dimensional chain of these 
unit cells as shown in Fig. [2j Such a billiard has finite 
horizon, i.e. it does not have unbounded collision-free 
trajectories. Note the aforementioned unit cell mirror 
symmetry in the transversal axis {x — 0). We focus on a 
few values of A2 in the range 0.3 < A2 < 0.60, where we 
found that the dynamics is strongly chaotic. Very small 
tori are observed for A2 — 0.38 and 0.45 in the longitudi- 
nal Poincare section over y = hi{x). For the others values 
of A2 considered the phase space looks fully ergodic to the 
naked eye. In all cases the connected ergodic component 
makes up at least about 95% of phase-space volume and 
the most important property for our purposes, exponen- 
tial decay of the instantaneous velocity auto-correlation 
{vx{t)vx{0)), was observed for all values of A2 explored 




FIG. 3: Velocity autocorrelation function C{t) — {vx{t)vx{0)} 
for an ensemble of 20000 particles in three billiards with A2 = 
0.6, 0.45, 0.3 (full, dashed and dotted lines respectively). The 
decay to noise level is quite fast in all cases, at around f ~ 15 
which is approximately 30 collisions. 



(Fig. [3]). This implies that the system exhibits normal 
diffusion, i.e., that 



X — lim 



x{t) - x{0) _ 



lim Xt 



(16) 



is a stationary Gaussian random variable, with zero mean 
and variance (x^) = Di the diffusion coefficient of the 
guide. In Fig. [4] we plot the distribution of the normal- 
ized displacement x for A2 = 0.3 showing remarkably 
good convergence to a normal distribution. This holds 
true for the other values of A2 that we consider. An- 
other way to support our conclusion that x is normal- 
distributed and that there is convergence for at least 
the first two moments 20 comes from the fact that 



i(litl)^ 



for large enough time as observed in Fig. 



[5| Note that, for instance, an anomalous diffusive system 
with ergodic phase-space but infinite horizon trajectories 
could have convergence in distribution to a Gaussian for 
an appropriate normalized displacement but would fail 
to satisfy this relation between first and second moments 
for any finite time approximation of x [21] . 



B. Bloch basis 

The propagating states of a periodic waveguide are 
given by the Bloch modes with unitary eigenvalues. The 
Bloch basis of a given waveguide can be defined by the 
solutions of a generalized eigenvalue problem as follows 
[22| . Consider only one unit cell between two infinite 
plane leads. The wave-function can be written as 

00 

i'{x,y)^Y.i^ti^) + c7{x))Mx,y) (17) 

where i E N labels the local transverse modes (j)i{x,y) 
which satisfy the waveguide boundary conditions at each 
X. In the particular case we are considering (hard-wall 
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by 



FIG. 4: Histogram of the normalized displacement Xt for a 
waveguide with yl2 = 0.3 at time t — 50000 using an ensemble 
of 10^ initial conditions. The best Gaussian fit is plotted over 
it. The inset shows the same histogram (dots) in log-scale 
where convergence to a Gaussian (full line) is more evident to 
be achieved even deep in the tails. 




FIG. 5: Finite time approximations of the diffusion coeffi- 
cient Di{t) = {xt) /t for the A2 = 0.3 cosine waveguide calcu- 
lated from the variance (second moment) of the sample's dis- 
placements (black line) and from the first moment using that 
f (l^tl)^ ~ i^t), as expected for Gaussians, (dashed line) for 
an ensemble of 10^ initial conditions. To obtain the diffusion 
coefficient Di we average D\{t) for times after the transient 
relaxation regime. The diffusion coefficient obtained in this 
way from the first and second moments of Xt are the same up 
to standard error. 



■ sm 



billiard chain) (f)„{x,y) ^ ^ 

h{x) — h2{x) — hi{x). Inserting Eq. (jlTl) in Eq. dl 



h(x) 




wave equation becomes a system of ordinary differential 
equations for the cf{x) which after suitable transforma- 
tion can efficiently determine transmission and reflection 
matrices [23ll2l]- Let (a~) be the right (left) going 
wave on the left lead and bf {b~) the right (left) going 
wave on the right lead (identified with cf{x) evaluated 
at the left xi and right Xr boundaries of the cell respec- 
tively) as illustrated in Fig. [61 If i and r (i' and r') are 
the left-to-right (right-to-left ^transmission and reflection 
matrices then outgoing waves are linked to ingoing waves 



-t'b- =a- 



(18) 



We can recast this (infinite) system of equations as 



Ml 



where 



Ml 



t 

-r I 



Mo 



b+ 

b- 



I -r' 
t' 



(19) 



(20) 



Now, we impose Bloch condition, namely that the wave- 
function is equal on both lead boundaries up to a complex 
factor A, 



(21) 

where A — exp{i9) — exp{iqL). Then, the Bloch basis 
is defined as the set of modes that satisfy (21), i.e. they 
are the solutions v„ of the unit cell's generalized transfer 
matrix eigenvalue problem 



Miv„ = A„M2V„ 



(22) 



The propagating Bloch modes are such that |A„| = 1; 
the other modes are evanescent and decay exponentially 
in the direction of propagation. In practice, this infinite 
system of equations must be truncated to a dimension at 
least equal to the number of propagating Fourier modes 
in the leads but in most cases it is necessary to include 
some evanescent Fourier modes to describe the field ad- 
equately. Note that the transmission and reflection ma- 



trices used in ( 20 ) are evaluated at the (right or left) en- 
try boundary of the unit cell and carry information from 
these decaying modes in the near field, which in general 
can couple between unit cells in the chain. Hence, these 
matrices do not satisfy the usual (energy) fiux conser- 
vation relations coming from the unitarity of the scat- 
tering matrix in the far field [25] . Moreover, including 
evanescent modes brings about infinitely fast decaying 
states (i.e. null eigenvalues A„) [55] making the matri- 
ces (20) non-invertible. In this case the transfer matrix 
M = M^^Mi cannot be defined. Only when the geom- 
etry of the unit cell is such that the coupling between 
evanescent modes is negligible it is possible to truncate 
t, r, t' and r' to the far field non-decaying modes and 
define the transfer matrix. 

Note that our method computes 6* as a function of k, 
i.e., we fix the energy E (or the wavenumber k) and then 
we obtain the set of Bloch quasi-momentum 0„ for which 
a propagative solution exists. With this method we can 
accurately compute the energy average at the left hand 
side of Eq.([7|. 
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FIG. 6: Illustration of a quasi-one-dimensional waveguide 
composed of one unit cell (gray) connected to two ideal 
leads going to infinity. The field (Eq. |17[ ) at the unit 
cell's boundaries is given by [cf{xo) c~(a;o)] = [af a^] and 
[cf{xi) c~(a;i)] — [bf . Solving the scattering problem 
in this geometry we can obtain the Bloch basis of the infinite 
periodic chain. 



IV. NUMERICAL RESULTS 

We now study numerically the semiclassical behavior of 
{NB{k)) for the cosine billiard chain. In order to compute 
the transmission and reflection matrices we employ the 
admittance multi- modal method described in |24| . where, 
for each value of k, the field is described using the local 
transverse modes (Fourier basis) of the waveguide trun- 
cated to a finite number of channels, as described in Sec. 



IIIB The number of propagating modes Nsik) is ob- 
tained by counting the solutions of Eq.(p2| with |A„| = 1 
and positive group velocity. As a function of fc, the above 
number fluctuates at scales Ak ~ 2tt/ (kAc) [IH]- We con- 
sider an average of the form 



{NB{k))r 



NB{k')fr{k - k')dk' 



(23) 



where fr{k) is a positive function of unit norm in dk 
with compact support in an interval of length rA/c with 
r 1. The dependence of {NB{k))r in the width of 
Jr[k) is briefly discussed at the end of this section but 
for what follows we drop explicit reference to r. 

In Fig. [t] we show {NB{k)) for k e [vr, ISOvr], for 
A2 = 0.3. For this A2 value, the phase-space does not 
show noticeable tori. In the k interval explored the 
number of open transverse modes in the Fourier basis 
ranges from 1 up to 150 (100 evanescent modes were 
used for all wavenumbers) . The dashed line shows a 
gaussian smoothed moving average with variance tt^ 
of NB{k) for a uniform sampling of wavenumbers with 
spacing Sk ^ 0.377r, which is bigger than the mean 
level spacing in the whole k range. The full line is 
the expected semiclassical result ( 15 ) adjusted with an 



additive constant. This constant appears because in the 
system under study {NB{k*)) = for some fc* > and 
we are not sufficiently deep in the semiclassical regime 
for its relative value to be negligible. NB{k) fluctuates 
very rapidly with k as expected but the agreement 
between the moving average and the semiclassically 
predicted curve is quite good in the whole wavenumber 
interval explored. We observe similar good results for 




FIG. 7: The dashed line is the gaussian smoothed moving 
average (with standard deviation ~ tt) of the immber of 
propagating Bloch modes Nsik) in the cosine billiard chain 
with A2 — 0.3 for a uniform k sampling with increments 
Sk = 0.371477r. The expected semiclassical result for {NB{k)), 
plotted in full line, shows good agreement with the moving av- 
erage. The four black points (with their associated standard 
error) were computed for the same geometry using a uniform 
sampling of 153 wavenumbers in an interval 300 mean level 
spacings wide, which constitutes a better statistics than the 
moving average (a sampling at least 15 times more dense). In 
the inset, the diffusion coefficient calculated classically (black 
dots) and the one obtained from the best fit following Eq. 
( |15[ ) for the quantum computation of Nsik) (circles) for five 
values of A2 is shown. In all cases, a good agreement with 
the expected semiclassical behaviour Eq. (151 is observed. 



other configurations of the unit cell with 0.6 > A2 > 0.3 
even when there are some small tori in the transverse 
Poincare sections. This is not surprising because the 
classical dynamics in those systems is also strongly 
mixing and for waves to resolve such small structures 
wavenumbers of order k ^ 20007r would be needed, a 
fairly large value compared to the range considered. 
In the inset of Fig. [7] we show the diffusion coefficient 
calculated classically (from Montecarlo simulations) and 
also the one obtained from the best fit following ( 15 ) for 



the quantum computation of NB{k) for five values of 
A2. 



As mentioned in the definition Eq. ( 23 ) of the mean 
number of propagating Bloch modes, the width of in- 
terval rAk over which we take the average is arbitrary 
and it is only assumed to be much larger that the mean 
level spacing Ak and classically small. In order to find 
out how relevant is this choice we calculate {NB{k))r for 
four different values of fc as a function of r G [2, 300] us- 
ing a uniform k sampling of 153 points around the four 
wavenumbers. In Fig. [8] we see that, after a transient 
phase in r '-^ 200, the value of the average reaches an 
approximately stable region in the four cases. The four 
black dots in Fig. [7] show this result for r = 300, also 
agreeing nicely with the expected semiclassical curve. 

We end this section remarking how close to the semi- 
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FIG. 8: {NB{k)}r as a function of r for four energies in the 
A2 = 0.3 cosine billiard. A uniform sampling of 153 wavenum- 
bers distributed in a 300 mean-level-spacing wide interval was 
used. In all cases after r ~ 200 the average fluctuations are 
notably reduced and the result for {NB{k)) is in agreement 
with the expected semiclassical result (see Fig. [7| indicat- 
ing that this is a good smoothing window to eliminate the 
spectrum fluctuations. 



classical regime the presented results are. In our calcu- 
lations X/W (the wavelength in units of the waveguide 
width) is as small as 0.012. We note that the calculations 
in [5j , which are also compared with semiclassical results 
for periodic billiards, were computed with the parameter 
\/W - 0.12. 



V. DISCUSSION AND CONCLUSIONS 

In section |ll] we closely followed Faure derivation of 
Eq.(13|. Here we provide an alternative argument al- 



lowing to understand better the connection with [51 [S] . 
Proceeding as in Sec. |llj we define the probability distri- 
bution of level velocity in natural units w 



p{w,E) 



so that 



27r ^ 



d0S{E - E„^g)6{w- 



dE. 



n.e ' 



Nb{E) = 2tt wp{w, E) dw. 



(24) 



(25) 



They showed this result for disordered systems and 
verified numerically its validity for chaotic systems in 
the semiclassical limit. It follows from their expression 
for k{uj,x) that the probability distribution of level 
velocity p{w,E) is a Gaussian with variance C{Q)/p\. 
Thus, from Eq.lps]), we obtain Nb{E) = ^2t:C{Q). 
The constant C(0) is system dependent and must be 
computed, for example, semiclassically. In our case 
it is simply related to the average of the square of 
the quantum velocity C(0) = ^f^(^^?i0) that is given 

leading to Nb{E) = \fl^^^- Thus, 



in Eq.(ll| 



Fame's result Eq.(13) can be seen as a consequence of 
universality in the parametric level correlations. The 
validity of the universal correlation function for periodic 
systems was studied in |B] by Dittrich et al. in the 
semiclassical limit. They computed K{Vt^ 0) by Fourier 
transforming their semiclassical expression for the form 
factor identifying a C(0) that agrees with the previous 
result. 

The universal correlation function k(uj,x) is different 
for GOE and GUE. For periodic systems, when x plays 
the role of the Bloch parameter, time reversal symmetry 
T is broken for all x besides three exceptional points 
[5]. Thus, when no further anti-unitary symmetries 
exist, fc(a;,x) moves between GOE and GUE at different 
parts of the band. In our case, the unit cell is invariant 
under the operator Sx that transforms x — —x (mirror 
reflection along the transversal direction). So, even 
though the unit-cell hamiltonian does not commute with 
Sx nor with T, it does commute with the operator TSx 
which is anti-unitary. Then, as x is varied, the system 
always belongs to GOE [28] . 



As a final comment, we note that Eq.(14) applies also 



to classical waves with dispersion relation uj — ck propa- 
gating in a wave-guide with speed c = 1. We can restore 
c by noticing that the distribution of group velocities 



Vn.e 



L 



de 



is related to our previous distribution by 



Pk{v) = \Pk{^)- With this new pk{v) the formula be- 
comes 



NB{k) 



2ncL 



vpkiv) dv 



(27) 



Then, as was noted in [27] , the distribution of level veloc- 
ities satisfies p{w, E) = liuVfp^a 4'K{w(j), (f)) with K{^1, (f) 
the autocorrelation function 



{p{E + n,e + cp)p{E,e))Efi 



Pe 



(26) 



of the energy density p{E,e) = J2n^i^ - En,e)- 
The average is over 9 and over a classically small 
energy interval around E and pe denotes this aver- 
age for p{E,9). According to Simons and Altshuler 
fc(w,x) = Kin,cl))/p- 

LO = psfl and X = V^^i^)^ with C(0) 



I; is a universal function where 
Pe{[- 



de 



with the normalization / dvpk{v) = / ^-Pfe(^) — 2ttAc 
for the group velocity distribution. The result for the 
number of propagating modes also applies for these 
classical waves. 

To conclude, we have studied the agreement between 
numerical computations and a semiclassical approxima- 
tion of the number of propagating modes in a diffusive 
periodic waveguide. We also showed the relation between 
{Nb{E)) and previously studied spectral properties of 
periodic systems. The semiclassical results used in our 
analysis is previous to the discovery of how to include 
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non-diagonal contributions to correlation functions. Nev- 
ertheless it was noted [HI [14] that their role in the com- 
putation of the fluctuations in velocity should be very 
limited. Indeed, besides the average justified by the rapid 
phases variation, the observable amplitude (the velocity 
in this case) fluctuates around its null average, further 
justifying the diagonal approximation. 
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